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ABSTRACT 



Aims. We construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegration 
of heavy elements such as Ru, which are believed to occur during superbursts in mixed H/He accreting systems. We use this network 
to simulate a detonation propagation, inside a mixture of C/Ru. 

Methods. As our reference, we use a full nuclear reaction network, including 14758 reactions on 1381 nuclides. Until the reduced and 
full networks converge to a good level of accuracy in the energy production rate, we iterate between the hydrodynamical simulation, 
with a given reduced network, and the readjustment of a new reduced network, on the basis of previously derived hydrodynamical 
profiles. 

Results. We obtain the thermodynamic state of the material after the passage of the detonation, and the final products of the combus- 
tion. Interestingly, we find that all reaction lengths can be resolved in the same simulation. This will enable C/Ru detonations to be 
more easily studied in future multi-dimensional simulations, than pure carbon ones. We underline the dependence of the combustion 
products on the initial mass fraction of Ru. In some cases, a large fraction of heavy nuclei, such as Mo, remains after the passage of 
the detonation front. In other cases, the ashes are principally composed of iron group elements. 

Key words. Hydrodynamics - Nuclear reactions, nucleosynthesis, abundances - Shock waves - Stars: neutron - X-rays: bursts 



1. Introduction 

Superbursts have been discovered by long-term monitoring of 
the X-ray sky, using instruments such as RXTE and BeppoSAX. 
Compared to normal type I X-ray bursts, superbursts are 1000 
times more energetics, with integrated burst-energies of approx- 
imately 10 42 ergs; they last 1000 times longer, from hours to half 
a day; and they reoccur every one to a few years. They are rare: 
only 13 events have been identified from 8 sources (for reviews 
see Kuulkers 2004, Cumming 120051 and references therein). 
They exhibit similarities with normal type I X-ray bursts, such 
as a rapid rise in the light curve, a quasi-exponential decay, and 
a hardening of the spectrum during the rise followed by a soften- 
ing during the decay. These attributes are well represented by a 
blackbody model, which has an effective temperature that grows 
during the rise, and decreases during the decay phase (Kuulkers 
2004). This suggests that superbursts, like normal type I bursts, 
are thermonuclear in origin (Cornelisse et al. 12000) . The current 
view is that superbursts are due to thermally unstable ignition 
of 12 C at densities of about 10 s - 10 9 g cm' 3 (Cumming I200T1 
Strohmayer & Brown [20021 Cummmg 120051 

Most superbursts have been observed in systems that accrete 
a mixture of H and He (Kuulkers 2004). In these H/He accretors, 
the 12 C mass fraction that remains after the combustion of H and 
He, via the rp-process (Wallace & Woosley 1981), is unknown 
(Schatz e t al. I200T1 |2003bl Woosley et al. [20041 . Cumming & 
Bildsten (1200 It showed that a 12 C mass fraction of 0.05 - 0.1 
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was sufficient for an unstable ignition; however, fits to the ob- 
served lightcurves suggest that the actual abundance is much 
larger (Cumming et al. 2006). 

Schatz et al. (2003) showed that, at the high temperatures 
reached in superbursts (T > 10 9 K), the photodisintegration of 
the heavy rp-process ashes influences the energetics of the det- 
onation. More details about this interesting mechanism are pro- 
vided below. 

To calculate energy production, a simplified nuclear reaction 
network is, in general, used in astrophysical hydrodynamics sim- 
ulations (e.g. Noel et al. 2007). In this study, we extend such a 
network to include a contribution by rp-process ashes. However, 
full nuclear reaction networks are intractable in hydrodynamical 
simulations, where some 10 7 calls to the nuclear reaction net- 
work, are required. This would be too time-consuming, in par- 
ticular in multi-dimensional simulations that will be considered 
in a future work. For this reason, we constructed a reduced nu- 
clear reaction network. This network and the detonation profiles 
are described in Sect. [2] Our results are discussed in Sect. [3] 



2. Nuclear reaction network extension 

In a mixed H/He accreting system, the detonation that leads to 
the superburst, is expected to propagate into the ashes of the rp- 
process, which occurred in the upper atmosphere. For simplic- 
ity, as in Cumming & Bildsten (2001), we represent the ashes 
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using a fiducial heavy nucleus, 96 RuQ We simulate the propaga- 
tion of a detonation, in a mixture of 12 C and 96 Ru, with mass 
fractions X»r u = 1 - X, 2 q, using the hydrodynamical algo- 
rithm described in Noel et al. (2007). This algorithm solves the 
adiabatic Euler's equations for compressible, non-viscous gas- 
dynamics with source terms. It is a finite-volume method in the 
spirit of the original MUSCL scheme of van Leer (1979). The 
algorithm is of second-order in the smooth part of the flow, and 
avoids dimensional splitting. Parallelization, which is required 
to solve computationally, the problem in hand, is based on the 
mpi library, as described in Deledicque & Papalexandris (2006). 
The equation of state accounts for partially degenerate and par- 
tially relativistic electrons and positrons. The ions are treated as 
a Maxwell-Boltzmann gas, and the radiation, considered to be 
at local thermodynamic equilibrium with the matter, follows the 
Planck law. Coulomb corrections to the equation of state are not 
included. In our simulations, these corrections would change the 
thermodynamic quantities of our ideal plasma by less than 2-3% 
(Fryxell et al. 2000), but will be included in future work. 

It remains untractable at the present time to perform the hy- 
drodynamical simulation with an extented network, so that a re- 
duced network needs to be constructed. For this purpose, we first 
extend the reduced network of thirteen species from 12 C to 56 Ni 
used in Noel et al. (2007), with the a-chain between 64 Ni and 
96 Ru. We adopt the nuclear data from the Brussels nuclear reac- 
tion rate library (BRUSLIB), based on published compilations 
of experimental reaction rates, and on the determination of reac- 
tion rates using the statistical Hauser-Feshbach model (Aikawa 
et al. 120051 Arnould & Goriely 2006). The photodisintegration 
rates (y, a) are calculated using the reciprocity theorem (Eq. (4) 
of Arnould & Goriely 2006). The resulting network is referred 
to as netO. 

Table 1. Initial conditions on the left- (column 2) and right-hand 
(column 3) sides of the discontinuity. 





Left 


Right 


P (gem J ) 


3.01 x 10 s 


10* 


r(K) 


4.46 x 10 9 


10 s 


v (cm s~') 


8.07 x 10 s 






Using netO, we simulate a detonation, propagating in a mixture 
X n Q = 0.2, X%r u = 0.8. The ignition conditions are the same as 
used for pure 12 C in Noel et al. d2007l > (see TableQ). The length 
of the domain is 1000 cm, the initial discontinuity is placed at 
x = 100 cm, where x is the distance from the left boundary of 
the domain, and the resolution is 1 cm (1000 numerical cells). 
We adopt a mixture X 5b ^ =0.1 and -^Ni = 0-^> at tne left 
of the discontinuity. The nuclear energy generation, tempera- 
ture, density and pressure profiles, at time t — 7 x 10 -7 s are 
presented in Fig. [1] in terms of the distance to the shock. The 
detonation velocity isZ)=1.19xl0 9 cm s _1 . Using the thermo- 
dynamic profiles presented in Fig.Q] we perform a full network 
calculation, including 14758 reactions on 1381 nuclides lying 
between the proton and neutron drip lines and with charge num- 
bers Z<50. Rates are taken from experiments, whenever avail- 
able, and otherwise from the BRUSLIB library. The derived nu- 

1 Choosing only % Ru implies naively that the electron captures in 
the high-density environment have not occured yet and that only % Ru 
is produced during the combustion phases preceding the superburst. 
Future detailled hydrodynamical simulation in X-ray burst may provide 
us with better initial conditions for the superburst. 
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Fig. 1. Nuclear energy generation (erg g _I s _1 ), temperature (K), 
density (g cirT 3 ) and pressure (erg cm~ 3 ) profiles for a detona- 
tion front in a mixture X n Q = 0.2, xY%r u = 0.8 at T = 10 8 K and 
p = 10 8 g cirT 3 . Z is the distance to the shock in cm. 




t[s] 



Fig. 2. Nuclear mass fractions of some species calculated with 
the full network including 14758 reactions on 1381 nuclides 
using the profiles of Fig. Q] The time in s is the time elapsed 
since the passage of the shock, assuming a constant velocity 
D = 1.19 10 9 cms" 1 . 



clear mass fractions, for some species are presented in Fig. [2] In 
Fig. [3] we compare the energy production rates of the reduced 
network netO, and of the full network. In this case the reduced 
network has reproduced the total energy production as well as 
the energy rate, of the full network. This favorable situation may 
not however, be possible, for detonations propagating in a mix- 
ture composed of a lower amount of 12 C. 
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Fig. 3. Network calculation of the energy production rate (erg 
g _I s _I ) and total energy (erg g -1 ) for a detonation in a mixture 
XnQ = 0.2, X,6R U = 0.8 using the profiles of fig.Q]for two nu- 
clear reaction network. Solid lines: full network including 14758 
reactions on 1381 nuclides; dotted lines: reduced network netO. 
The time in s is measured since the passage of the shock, assum- 
ing a constant velocity D = 1.19 10 9 cm s _1 , 
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Fig. 4. Same as Fig.Q]but for a mixture X u q — 0.1, X»r u =0.9. 



In a second step, we simulate a detonation propagating in a mix- 
ture XnQ =0.1, ^»r u = 0.9. The ignition conditions are similar 
to those used above. The nuclear energy generation, tempera- 
ture, density, and pressure profiles, at time t — 7 x 10~ 7 s are 
presented in Fig. [4] in terms of the distance to the shock. The 
detonation velocity isD=1.16xl0 9 cm s _I . Using the thermo- 
dynamic profiles of Fig. |4] a full-reaction network calculation is 
performed. The derived nuclear mass fractions, for some species, 
are presented in Fig. [5] 

We compare the energy production rates, of the reduced net- 
work netO, and of the full network (Pig. |6j>. We see that this case 




i(v* io a irr irv 6 io 5 
t[s] 

Fig. 5. Same as Fig.[2]but for a mixture X 12 q =0.1, ^>6R U = 0.9, 
using the profiles of Fig.|4]and assuming a constant velocity D = 
1.16 10 9 cms" 1 . 



is characterized by an endothermic phase, due to the enhanced 
initial abundance of 96 Ru, and the significant contribution of the 
photodesintegration reactions, to the total energy balance, at an 
early stage. This phase of the energy generation profile obtained 
using the full network, is not reproduced by the reduced network 
netO. We therefore construct a more sophisticated reduced net- 
work, called netl. Let us describe the methodology used in its 
construction. In the reduced networks, only the capture or emis- 
sion of a-particles are taken into account. However, neutrons and 
protons can be produced significantly by photo-reactions, and be 
recaptured. In particular, the nucleus (Z-2.N-2), where Z is the 
charge number, and N the number of neutrons, can be produced 
by the photodesintegration of the nucleus (Z,N), following six 
additional paths, complementing the (y, a) one; these include 
the emission of two neutrons and two protons. The characteristic 
timescale, for each of the seven possible paths is given by 

i-i _ i-i 

_ A (y,a)> 

^ = C/ z > + C)< z ' N ~ l ) + K.p^ N - 2 ) 

+^)(Z-l.tf-2), 
A 2 1 = A^ n) (Z, N) + A^ p) (Z, N-l) + A^ n) (Z - 1, N - 1) 
+A^ p) (Z-l,N-2), 
= *(y\ n) (Z,N) + A- y \ p) (Z,N - 1) + A-^JZ -\,N-\) 

+X£jZ-2,N-l\ 
= ^, P) (Z, AO + A^JZ - 1, AO + A- y \ n) (Z -l,N-l) 
+A^ p) (Z-l,N-2), 
A- 5 l = A^ p) (Z, N) + A^ n) (Z - 1 , N) + 17 y \ p) (Z -l,N-l) 

+Al y \ n) (Z-2,N-l), (1) 

so that the total photo-reaction rate of (Z,N), leading to (Z-2,N- 
2), can be approximated by the effective rate 

Ae f f = A + Ai + A 2 + A 3 + A 4 + A 5 + A 6 . (2) 
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The reverse a-capture rate can then be calculated using the de- 
tailed balance equation (Eq. (4) of Arnould & Goriely 2006). 
While /lo, describing or-photoemission, is the only path included 
in the netO approximation, the neutron and/or proton photoe- 
mission could potentially dominate in the case of light (Z<40) 
species, depending on the nuclei involved. In the neutron- 
deficient region (especially along the neutron-shell closures), the 
proton-photoemission is dominant, while close to the stability 
line the neutron-photoemission usually dominates. Any of these 
seven paths, however, can be impeded to a more or less large 
extent by the reverse nucleon captures that would modify the 
relative significance of a given path. This modification depends 
on the thermodynamic and composition conditions, as well as 
on the nuclear cross sections. To determine which paths are im- 
portant, we use the following method. We first calculate the det- 
onation profiles using netO. On the basis of the temperature and 
density profiles (for instance those of Fig. [4] for X I2 q = 0.1, 
X^Ru = 0.9), the full reaction network is solved. This calcula- 
tion allows to determine the dominant paths and to construct a 
reduced network with new effective reaction rates but the same 
reduced number of nuclei which is introduced in the hydrody- 
namical algorithm. We perform the detonation simulation using 
this new reduced network, and the resulting thermodynamical 
profiles are used again in a full network calculation. We iterate 
between a hydrodynamical simulation with a given reduced net- 
work, and the readjustment of a new reduced network, on the 
basis of the previously-derived hydrodynamical profiles, until 
the reduced and full networks converge to achieve high accu- 
racy in the energy production rate. Note that no strict mathe- 
matical criteria have been used to estimate the quality of the 
newly-developed reaction network. In contrast, we estimate that 
the new network must fulfill two major conditions: it must re- 
produce the right order of magnitude (roughly within a factor 
of two) of the total amount of energy produced, and the possi- 
ble appearance of endo thermic phases. Indeed any endo thermic 
phase can play a role of first importance on the detonation prop- 
erties. In particular pathological detonations may occur (Fickett 
& Davis [T979l Khoklov [19^1 Sharpe [19991 . 

In contrast to the reduced network used in Noel et al. (12007b . 
which provided an energy source typical of explosive helium and 
carbon burning in the absence of hydrogen (Gamezo et al. 1999), 
our reduced network depends on the detonation profiles, and ini- 
tial composition. For each new detonation simulation, we con- 
struct the adapted reduced network using the iteration procedure 
described above. 

Convergence is achieved in the case of X 12 q =0.1, = 
0.9, with the following reduced network, called netl, in which 
the equivalent reaction rate, adopted for each photodisintegration 
from 92 Mo to 64 Ni, and from 56 Ni to 16 Q is 



Aq + Ai + A(„ 



(3) 



and for Ru, A ne ,\ — Aq. The inverse a-capture rates are 
calculated using Eq. (4) of Arnould & Goriely (2006), and 
the five remaining rates for I2 C( 12 C,a) 20 Ne, 12 C( 16 0,a) 24 Mg, 
16 0( 16 0,a) 28 Si, 4 He(2a,y) I2 C, and 12 C(y, 2a) 4 He, are similar 
to those used in netO. 

Figure [6] compares the energy production calculated using 
netO, netl, and the full network. Using the profiles obtained with 
netO (Fig. @J, we perform a full network calculation, and com- 
pute the energy production rate. The same profiles are used for 
the reduced network netO, and the energy production rates are 
compared. We see that the endothermic part of the profile is not 
reproduced by the reduced network netO. A revised reduced net- 
work is therefore developed, netl, which simulates more accu- 
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Fig. 6. Same as Fig.[3]but for a mixture X 12 q =0.1, X>6R U = 0.9, 
assuming a constant velocity D = 1.16 1 9 cm s~'. and for 
three nuclear reaction network. Solid lines: full network includ- 
ing 14758 reactions on 1381 nuclides; dotted lines: reduced net- 
work netO; dashed lines: reduced network netl. 
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Fig. 7. Same as Fig. |4] but with netl. 



rately, the full network energy production rate. However, new 
detonation profiles, using the hydrodynamical algorithm have to 
be calculated with this new network. They are shown in Fig. [7] 
at time tf - 1.6 x 10~ 6 s. The domain length is of 2000 cm, and 
the resolution is of 1 cm. The initial discontinuity is positioned 
at x = 200 cm. The detonation velocity is D = 1 .09 x 10 9 cm s~' . 

Using these new thermodynamic profiles (Fig. |7J, network 
calculations of the energy production are made for the full net- 
work, and the reduced network netl, which are compared in Fig. 
[8] One sees that the reduced network reproduces satisfactorily 
the energy production given by the full network. We can there- 
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Fig. 8. Same as Fig. [6] but with the profiles of Fig. [TJand D 
1.09 10 9 cms- 1 . 




Fig. 9. Same as Fig. [5] but with the profiles of Fig. [7j and D 
1.09 10 9 cms- 1 . 



fore consider netl as a suitable reduced network for the simula- 
tion of a detonation for a mixture X xl c =0.1, ^%r u = 0.9, and 
the initial conditions shown in Table [1] Adopting the thermody- 
namic profiles shown in Fig. [7] a final calculation is performed 
using the full network. The derived nuclear mass fractions for 
some species, are presented in Fig. [9] A comparaison between 
Fig. [5] and Fig. [9] reveals the sensitivity of the final composi- 
tion to the choice of nuclear reaction network. For instance, netO 
produces much less 90 Zr and 92 Mo than netl, at the end of the 
reaction zone. 

Far more than two reduced networks were tested in the iter- 
ation procedure. Figure [10] displays the total energy obtained, 
using five different reduced networks, for the thermodynamic 
profiles of Fig. [4] It is clear that netl reproduces most closely 
the global energy obtained by the full network simulation. 

Note that netl includes photoneutron and photoproton paths 
that are found to reproduce globally all possible paths, when in- 
cluding neutrons and protons in the network. Both paths (path 
1 and path 6 in Eq. [TJ are required to simulate the bottle- 
necks which are related to the neutron (N=20, 28, 50), or pro- 
ton (Z=20, 28) magic numbers, which are crossed by the pho- 



Fig. 10. Same as Fig. |6]but for five different networks and only 
the total energy produced is presented. The network net2 is iden- 
tical to netO for all species between 4 He and 56 Ni, and uses A e ff 
for the rates of all reactions between 64 Ni and 96 Ru. The net- 
work net3 uses A e ff for the rates of all photo-reactions between 
16 and 96 Ru, and the netO rates for the other reactions. 



tonuclear flow. All nuclei included in the reduced network are 
even-even, such that if a first (y, ri) reaction is possible, the re- 
sultant (Z,N-1) nucleus will certainly photo-emit a neutron. The 
same holds for proton emission. The equivalent rate A net \ in Eq. 
[3] simulates the possible inclusion of neutrons and protons in the 
network. The inverse effective c-captures rate, that we deduced 
from detailed balance expressions, similarly enable to simulate 
the possible neutron and proton captures, without having to in- 
clude neutrons and protons explicitely in the network. Finally 
note that for 96 Ru, we do not take the (y,n) and (y,p) contribu- 
tions from Eq.[3j into account to counterbalance the neutron and 
proton captures by 96 Ru, at early time when the matter is almost 
entirely made of 96 Ru. 

We emphasize that the reduced network considered for the 
detonation of the mixture of X I2 q = 0.1, ^%r u = 0.9 cannot 
be generalized automatically to cases characterized by differ- 
ent initial conditions. The construction of the reduced network 
does not only depend on the temperature and density time evolu- 
tion, but also on the specific initial composition. Different initial 
mass fractions of the existing species, or a different elemental or 
isotopic composition, would require the construction of a new 
network, with new optimum effective reaction rates. As an ex- 
ample, we present a second case study, in which X I2 q = 0.05, 
X^Ru = 0.95, as considered by Cumming & Bildsten (2001). 
The detonation profiles, in this case, are shown in Fig.[TT] using 
netl. The domain length is 1000 cm, the resolution is 1 cm, the 
initial discontinuity is placed at x = 100 cm, and the detonation 
velocity is D = 9.88 X 10 8 cm s" 1 . We see in Fig.[12]that netl 
does not reproduce the endothermic phase. A more suited net- 
work needs to be defined, for example, by considering a linear 
superposition of photorates, as in Eq.|2] that is different for each 
nucleus included in the network. Wheter such a scheme should 
be developed depends on the outcome of realistic modelling, of 
the initial composition of superbursts, in heavy elements. 
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Fig. 11. Same as Fig.@but with X llQ = 0.05 and X %Ru = 0.95. 
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Fig. 12. Same as Fig. [6] but with the profiles of Fig. [TT] D = 
9.88 10 8 cm s- 1 and X uC = 0.05, X 96Ru = 0.95. 



3. Conclusion 

We have shown that a simple extension of an c-chain, might 
not be sufficiently accurate to reproduce the energy production 
rate of a detonation in a mixture that contains too much heavy 
elements. We have developed a new methodology to construct 
a reduced network, with adapted effective reaction rates, based 
on an iteration procedure between hydrodynamical simulations, 
with a chosen reduced network, and a comparison with full net- 
work calculations. In particular, we have constructed a new re- 
duced network that reproduces satisfactorily, the energy gener- 
ation rate of a reference network, including all reactions on all 
nuclides, probably involved in the specific case of the propaga- 
tion of a detonation at T — 10 8 K, and p = 10 8 g crrT 3 , in a 
mixture X 12 q =0.1, = 0.9. Adopting such a reduced net- 



work, we have simulated detonation profiles. These profiles are 
used in a full network simulation to compute the associated nu- 
cleosynthesis. One sees that the presence of 96 Ru decreases the 
total energy of the detonation, due to an endothermic phase at 
early stages. For this reason, the propagation velocity of the det- 
onation decreases. For X I2 q =0.1, ^9&r u = 0.9, essentially only 
iron-group elements remain at the end of the reaction zone, with 
some addition of N=50 elements, such as 90 Zr. This fact can in- 
fluence the neutron star crust composition and properties, such 
as its conductivity, or its neutrino emissivity, as well as the igni- 
tion conditions for the superbursts (Weinberg & Bildsten 2007 
Gupta et al. 120071 Haensel & Zdunik [l990i l. 

Interestingly, it can be seen that for > 0.9, the deto- 

nation goes through a short endothermic phase (see Figs. l7lfT2l 
that cannot be described by netO (Fig. |4j. The photodesintegra- 
tions by a emission are endothermic reactions, and only the re- 
combination of a particles up to Ni compensates for the energy 
loss. The energy is therefore released by the captures or recom- 
binations of the protons, neutrons and or-particles, liberated by 
the photodesintegrations. Furthermore a comparaison between 
Fig.|4]and Fig.|7]reveals that the temperature reached during the 
detonation is lower using netl, than netO, which influences the 
nucleosynthesis. 

Computationally mixed C/Ru detonations are relatively easy 
to handle since the length-scales on which the different elements 
burn are comparable. This contrasts with the case of pure car- 
bon detonations where a large variety of length-scales is ob- 
tained (Noel et al. 2007 ), and where multiple simulations at dif- 
ferent resolutions are required to fully resolve the detonation. 
In the C/Ru case, carbon and ruthenium burn on « 10 cm, and 
the total reaction length is as 1000 cm. The length-scales span 
only three orders of magnitude compared to six for pure car- 
bon detonations. This is an encouraging feature for future multi- 
dimensional calculations. 
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